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ABSTRACT 

In this paper we investigate the detectability of a habitable-zone exomoon around 
various configurations of exoplanetary systems with the Kepler Mission or photom- 
etry of approximately equal quality. We calculate both the predicted transit timing 
signal amplitudes and the estimated uncertainty on such measurements in order to 
calculate the confidence in detecting such bodies across a broad spectrum of orbital 
arrangements. The effects of stellar variability, instrument noise and photon noise are 
all accounted for in the analysis. We validate our methodology by simulating synthetic 
lightcurves and performing a Monte Carlo analysis for several cases of interest. 

We find that habitable-zone exomoons down to 0.2M(^ may be detected and ^ 
25, 000 stars could be surveyed for habitable-zone exomoons within Kepler^s field-of- 
view. A Galactic Plane survey with Kepler-class photometry could potentially survey 
over one million stars for habitable-zone exomoons. In conclusion, we propose that 
habitable exomoons will be detectable should they exist in the local part of the galaxy. 

Key words: techniques: photometric — planets and satellites: general — planetary 
systems — occultations — methods: analytical 



1 INTRODUCTION 

With the successful launch of Kepler, a mission specifi- 
cally designed to detect habitable Earth-like planets, a new 
age in extrasolar planetary science is about to commence. 
Kepler is a mission designed to detect the transit of an 
Earth across the Sun with its highly s ensitive pho t omet - 
ric c amera; more detai ls can be found in iBasri et al.l (|2005| ) 
and iKoch et all (|2007l). as well as on the mission website 



(http://www.kepler.nasa.gov/sci 



In this paper, we evaluate the range of exomoons that 
the Kepler Mission or _R'ep/er-class photometry (KCP) could 
detect through transit timing effects, with particular at- 
tention to habitable-zone exomoons. We find that habit- 
able exomoons down to Q.2M^ should be detectable with 
the expected performance of Kepler. The theory behind de- 
te cting exomoons through timin g effects was first proposed 
bv lSartoretti fc Schneider l|l999ll and has been r ecently ad- 
vanced bv lKippind l|2009ah and lKippin3 (l2009bl ) and these 
models will be utilised throughout our analysis. 

The underlying principle is that an exomoon should 
induce transit-time variations (TTV) and transit-duration 
variations (TDV) on the host planet. These effects are pre- 
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dieted to exhibit a 7r/2 phase difference providing a unique 
exomoon signature. In this work, we consider that both ef- 
fects must be detected in order to claim a moon is present. 
Additionally, measuring both timing effects allows for the 
determination of the exomoon's mass and orbital distance. 

The question as to whether current telescopes can de- 
tect exomoons could have far-reaching implications for the 
future of exoplanet science. Potentially, exomoons could be 
common habitable environments in the galaxy and would 
therefore be of great interest to astrobiologists. In this paper, 
we find that exomoons are indeed already quite detectable 
with KCP. We emphasise the use of /("epZer-class photome- 
try due to the increasingly impressive results being obtained 
from the groun d which are matching space-based photome- 
try, for examp le i Johnson et al.l (|2009l ). Furthermore, ground- 
based observations are often more ideally suited for transit- 
timing studies due to the fewer constraints placed on the 
system, such as telemetry-limited data-download speeds. 



TRANSIT TIMING EFFECTS DUE TO AN 
EXOMOON 



ISartoretti fc Schneider! (| 19991 ) predicted that an exomoon 
should induce transit time variations (TTV) on a host exo- 
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planet by virtue of the fact that the planet and moon orbit 
a common centre of gravity. More rece ntly, the t heory of ex- 
omoon detection has been refined in iKippind (|2009ah and 
iKippingj (|2009bl ). According to these papers, there should 
exist two principal transit timing effects: 

(i) Transit time variation (TTV) caused by the position of 
the planet oscillating around the common centre-of-gravity 
of the planet-moon system. 

(ii) Transit duration variation (TDV) caused by the ap- 
parent velocity of the planet increasing and decreasing as 
it moves around the centre-of-gravity of the planet-moon 
system. 

TTV and TDV must be detected in order to differ- 
entiate between a signal due to an exomoon or, fo r ex- 
ample, a perturbing pla net (see lAgol et al] l|2005h and 
iHolman fc Murravl l|2005h '). This is possible since TTV and 
TDV are predicted to be 7r/2 out-of- phase for an exomoon, 
thereby producing a unique detection signature. 

We briefly mention that other techniques for detect- 
ing exomo ons h ave b een proposed in the form of mi- 
crolensing, iHanl (|2008l ). lightc urve distortions, ISzabo et al.) 
l|2006l '). planet -moon eclipses, Cabrera fc Schneidej I 2007^ ■ 
pulsar timing, iLewis et al.l (|2008D and very weak distortions 
in the Ross i ter-M cLaughlin effect of a transiting planet, 
ISimon et al] (|2009| ). 



3 MODELLING THE DETECTABILITY OF 
EXOMOONS 

3.1 Confidence of detection 

In order to explore a large range of parameter space, it is 
more convenient and efficient to employ analytic expressions 
rather than repeated individual simulations for thousands 
of different scenarios. In order to proceed we need general 
expressions for the following: 

(i) TTV & TDV signal amplitudes 

(ii) Transit mid-time and transit-duration errors 

(iii) Confidence of detection, based on signal-to-noise 

The TTV and TDV root mean square (r.m.s.) ampli- 
tu des may be easily calc ulated us i ng the e quations presented 
in iKippin j ([2009a) and iKippinS (|2009bh . However, we also 
require expressions for the timing errors, which are critical 
in evaluating the signal-to-noise. 

To address this, we will use the analytic expressions for 
the uncertainty on the mid-transit time {tc) and duration 
(T) as derived bv lCarter et al.l (|2008l ) using a Fisher-analysis 
of a trapezoid-approximated circular-orbit lightcurve. These 
expressions consider the transit duration to be defined as the 
time taken for the pl anet to move between contact points 1.5 
to 3.5 (see figure 1 of ,Cartero t al. (2 0081)). This is in contrast 
to the definition used bv [Seager fc Mallen-OrnelasI l|2003l ). 
for example, who consider It as being the duration between 
contact points 1 & 4 and as t he duration between con- 
tact p oints 2 & 3 (see figure 1 of ISeager fc Mallen-OrnelasI 

diooi)). 

For the purposes of TDV measurements, the primary 
requirement is to use a measure of transit duration which 
has the lowest possible un certainty. By calculating the co- 
variances of the lightcurve, ICarter et al.l l|2008f ) were able to 



show that T can be calculated more precisely than either 
tr or tp and so we select T as a robust duration parame- 
ter to explore the TDV effect. The uncertainties on transit 
depth, d, t ransit duration, T, a nd mid-transit time, tc, were 
derived bv lCarter erall (|2008l ') to be: 



w^ = dyivr 



(1) 

(2) 
(3) 
(4) 



where r is the ingress/egress duration, Tph is the photon 
collection rate, and d is the transit deptlp. 

These expressions do not hold for a poorly sam- 
pled ingress or egress and therefore we assume a ca- 
dence of 1 minute, corresponding to Kepler's transit- 
timing/ asteroseismol ogy mode. 

The equations of lCarter et al.l (|2008h require the ingress 
duratior0, r, and the transit duration, T, as inputs. In or- 
der to compute these val u es we will use the expressions of 
ISeager fc Mallen-O rneia^ lj2003l ') modified for the time be- 
tween contact points 1.5 fc 3.5. Note, that we do not requir e 
a more elaborate model, such as that of IKippingj l|2008h , 
since we are dealing with circular orbits only. For the Kepler 
Missio n or KCP, we emplo y the sam e estimate for 
that of lBorucki et al. 1 l|2005l ) (B05) and lYee fc Gaudil (|2008h : 



- ph 



6.3X 10«hr-M0-°-"('"-i'\ 



(5) 



where m is the apparent magnitude. 

For a normal transit depth observed n times, the con- 
fidence, C, to which the transit is detected, in terms of the 
number of standard deviations, is defined by B05 as: 



C(photometric) = —y/n. 



(6) 



The transit timing signals due to an exomoon are peri- 
odic in nature and so require a different detection method. 
Typically, this problem is tackled by searching for significant 
peaks in a per iodogram, as often e mployed for radial velocity 
searches (e.g. iButler et af] (|2002l )). However, this approach 
is less useful for transit timing effects due to an exomoon 
since the frequency we are trying to detect will always be 
m uch higher than the sampling frequency, as pointed out 
bv lKippind (|2009al ). Since we are always below the Nyquist 
sampling rate, then only a range of possible harmonic fre- 
quencies are retrievable and a hence a set of possible exo- 



moon masses. 



We propose here that the ideal method is to search for 
statistically significant excess variance and then use the x^- 
distribution to calculate the confidence of signal detection. 
In order for this method to be applicable, we require a) that 
the uncertainty estimates are robust and accurate; and b). 



^ We choose to cliange the notation shghtly from the original 
paper to avoid confusion with other parameters we use here 
^ Defined as the duration between contact points 1 & 2, which is 
of course equivalent to the egress duration for a circular orbit 
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that the period of the signal may be derived from ampli- 
tude information alone (which may then be compared to a 
periodogram to further refine the frequency). 

The first of these requirements can be seen to be val id as 
several investigations have verified. iHolman et al.l (|2006l ) de- 
rived the uncertainties of the mid-transit times for four tran- 
sits of XO-lb using three different methods: i) Ax^ = 1 per- 
turbation of the best-fit; ii) Monte Carlo bootstrapping; iii) 
Markov-Chain Monte Carlo (MCMC). The authors found 
that all three methods produced very similar uncertainties, 
which implies the uncertainty estimate s are highly robust . 
Another example we point out is that of lCarter et al.l ()2008l ) 
who showed that the uncertainties derived using an MCMC- 
analysis were very similar to those predicted using analytic 
arguments. 

Th e second of our requirements is validated bv lKippind 
l|2009al ). where it was shown that the ratio of the TTV and 
TDV signal amplitudes may be used to obtain the period of 
the exomoon. This period may then be compared to the set 
of possible harmonic frequencies derived from a periodogram 
in order to obtain a highly reliable estimate. We therefore 
conclude that a search for excess variance is the most appro- 
priate strategy for searching for exomoons through transit 
timing effects. The confidence of detection may be found 
by integrating the probability density function of the x^- 
distribution. 



■ n = M / Pp where M is the mission duration and Pp is 
the period of transiting planet. 

■ At least three transits are needed to detect both a 
planet and a moon. 

In most of the cases we will consider, many more than 
three transits will be detected and three can be seen to be 
the limiting case for GOV stars, where the habitable zone is 
sufficiently distant to only permit three transits in a 4-year 
timespan. Although statistically speaking three transits is 
sufficient, there is a risk of an outlier producing a false posi- 
tive. We therefore consider detections of habitable exomoons 
in early G-type star systems to be described as 'tentative', 
whereas once four transits are detected, for stars of spec- 
tral type G5V and later, this risk can be considered to be 
reduced. 

The nominal mission length of Kepler is 3.5 years and it 
may be extended to up to 6 years, which justifies our choice 
of 4 years of transit timing observations. A ground-based 
search achieving KCP may easily be operational for 4 years 
or more. We choose 8-(t as the signal detection threshold 
since this is the same as that used by Kepler. The second 
signal may be detected to lower significance since it is only 
used to confirm the phase difference between the two and 
also derive the exomoon period. 



,r Z'"" T*"/^)"^ exp""'''^ 1 
C(timing)^V2erf--[l-/^^ 2^/^r(n/2) H' 

where n is the number of transits observed and is the 
observed value of , given by: 

"'="(l+A2)' (8) 

whe re 5 is the r.m.s . amplitude of the transit timing signal 
(see lKippinj l|2009al )') and A is the uncertainty on the mid- 
transit time/transit duration. 

Integrating and making the above substitution we have: 

C7(timing) = Vlerf-^ [i _ qg, ^ (l + ^) }] , (g) 

where Q{a, 6} is the incomplete upper regularized Gamma 
function. We summarise our assumptions below: 

■ Only one exomoon exists around the gas giant exo- 
planet of interest. 

■ The moon and planet are both on circular orbits and 
the moon's orbit is prograde. 

■ The moon's orbital plane is coaligned to that of the 
planet-star plane which is itself perpendicular to the line-of- 
sight of the observer, i.e. i = 90°. 

■ If a planet is within the habitable zone, then any moon 
around that planet may also be considered to be habitable. 

■ A transiting planet must be detected to 8-cr confidence 
to be accepted as genuine. 

■ An exomoon must be detected through either a) TTV 
to 8-cr and TDV to S-a confidence or b) TTV to 3-cr and 
TDV to 8-(T confidence, in order to be accepted as genuine. 

■ The Kepler Mission or KCP will be used in high ca- 
dence mode for the transit timing of a target of interest for 
~ 4 years. 



3.2 The total noise 

The expressions of ICarter et al.l (|2008l ) only consider shot 
noise through the Vph parameter. However, if we assume 
that the impact on ctt, ot^ and Od are approximately equiv- 
alent for additional uncorrelated noise and for correlated 
noise, then we may simply modify W to absorb the effects 
of red noise. For uncorrelated noise we may add the addi- 
tional sources of noise in quadrature. Note, that this is the 
same treatment utilised in the design technical documents 
for Kepler, for example see BOS. 

In general, there are expected to be three major types 
of noise present in the Kepler data in the form of shot noise, 
instrument noise and stellar variability. Instrument noise is 
due to a vari ety of effects and h as been modelled in depth 
by B05 (and ' Koch et al] l|2004h ') to quantify its effect as 
a function of magnitude. With all three noise sources, we 
modify W to W' , given by: 

T^4/5^+^' (1°) 

where / is the instrument noise and 5* is the stellar variabil- 
ity. 

7 is a function of magnitude which may be calculated 
using the model of BOS and we show all three noise sources 
plotted as a function of magnitude in figure 1. We assume 
a constant value for stellar variability of 10 ppm across 
all spectral types, a reasonable assumption, given that 6S- 
70% of F7-K9 main-sequence stars in the Kepler field are 
likel y to have similar or l ower intrinsic variability than the 
Sun (|Batalha et all (|2002h : BOS) on timescales important to 
transit detections. We also note that this equation is equiva- 
lent to the formulation used in the original technical design 
papers for Kepler, for example see equation (1) of BOS. 
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9 10 11 12 13 14 15 16 

Apparent magnitude 

Figure 1. Noise sources predicted to affect Kepler photometry 
as a function of magnitude. Instrument noise is dashed, photon 
noise is dot-dashed, stellar variability is dotted and the total is 
solid. Values obtained from Bill Borucki in personal communica- 
tion. 



3.3 The habitable zone 

We choose to consider a moon-hosting gas-giant exoplanet 
around a variety of main-sequence stars as shown in ta- 
ble 1. For each star we calculate the habitable zone or- 
bital distance, auab, to be defined as the distance where a 
planet would receive the same insolation as the Earth. This 
straightforwardly permits a reasonable estimate of the hab- 
itable zone for each star type. For a more in-dcpth consider- 
ation of habitability of exomoons see Williams et al.. (,1997iV 



ab = V^./LqAU. 



(11) 



is the same assumption used for Kepler^s ability to detect 
Earth-like planets since very active stars will be too vari- 
able for the detection of such bodies. In table 1 we list the 
different star properties. 



3.5 Properties of exomoons 

Although no exomoons have yet been discovered, it is pos- 
sible to calculate the range of ex omoons which are dynam - 
ically stable around each planet. iBarnes fc O'BrienI (|2002| ) 
addressed this problem and developed a set of analytic ex- 
pressions, which can be shown to provide excellent agree- 
ment to numerical simulations, for the stability of exomoons 
around exoplanets. Assuming an exomoon has to be stable 
for at least 5 Gyr, we are able to calculate the maximum al- 
lowed ex omoon mass in each case example (see equations (7) 
to (9) of lBarnes fc O'BrienI |20o5)). We assume Jupiter-like 
values for the tidal dissipation factor, Qp — 10^, and for the 
tidal L ove number, k^-p = 0.51, as used bv lBarnes fc O'BrienI 
(|2002l ). 

We also are able to estimate the range of allowed values 
for the planet -moon separation, in u nits of Hill radii, which 
we label as ^. jPomingos et al.l (|2006| ) presented the relevant 
expressions, which again can be shown to provide excellent 
agreement to numerical simulations. Using their equation 
(5), we are able to estimate the maximum distance at which 
an exomoon is stable for prograde orbits. 

For the minimum distance, we calculate the Roche limit 
of the planet in all cases. If this value is greater than 2Rp 
we use the Roche limit as the minimum distance, otherwise 
2Rp is adopted. We assume that there is no reason for a 
moon to exist at any particular value of ^ and thus the prior 
distribution is flat. 



For each planet-moon system we consider, the period 
of the transiting planet is calculated using Kepler's Third 
Law: 



Phab = 27r 



G{M, + Mp + Ms) ' 



(12) 



4 LIGHTCURVE SIMULATIONS 
4.1 Lightcurve generation 

A criti cal assumption in t his paper is the use of the equa- 
tions of lCarter et al.l (|2008l ) for the uncertainties on the tran- 
sit duration and mid-transit time. The authors tested their 
expressions using synthetic lightcurves and an MCMC fit- 
ting procedure as well as a numerical Fisher-analysis. They 
report excellent agreement between their expressions and 
the derived errors but find the greatest departure for poorly 
sampled ingresses and near-grazing transits. Our choice of 
assumptions avoids both of these issues. 

We note that these tests were run for shot noise only 
and in this paper we consider the effects of both instrument 
and stellar noise. Having modified W to W' to account for 
these additional noise sources, we choose to test this modi- 
fied formulation through generation of synthetic light c urves . 

In this work, we will use the iMandel fc Agoll (|2002l ) 
code to generate lightcurves acco unting for lim b darkening 
based upon a non- linear law from lClaret] (|2004l '). We use the 
ISeager fc Mallen-Ornelas (,200l) expressions to compute the 
motion of the planet since we are dealing with circular orbits 
only. We generate the lightcurves for a Neptune, Saturn and 
Jupiter-like planet in the habitable zone of a rriKcp = 12, 
Guidelines available from |http://keplergo.arc.nasa.gov/proposal.html |G2V Sun-like star with i — 90° and e = 0. Each lightcurve 



We choose to work in the time domain, rather than the 
orbital-distance formulation, since a major limiting factor in 
our study is the Kepler Mission duration. 



3.4 Properties of host star 

In our analysis, we will only consider single main-sequence 
stars which offer the best potential for hosting habitable en- 
vironments. We will consider spectral types from M5V to 
FOV and assume for each an ap proximate m ass, radius, and 
effective temperature as given bv lCoxl (|2000l ). and a luminos- 
ity derived from data therein. We use the Kepler bandpass 
to calculate the absolute magnitude of these star^. 

For each stellar type, we assume the stars are not young 
and may be considered to be slow rotators. Sin ce stellar vari- 
ability is correlated to rotational period (see jPorren et al.l 
l| 19941 )1. we therefore limit ourselves to quiet stars. This 
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Table 1. Properties of stars used in our calculations. Values taken from I Cod \200A) . Absolute magnitudes in the Kepler bandpass 
calculated using guidelines on the mission website. 



Star type 


M./M0 


R,/Rq 






Mkcp 


M5V 


0.21 


0.27 


0.0066 


3170 


11.84 


M2V 


0.40 


0.50 


0.0345 


3520 


9.49 


MOV 


0.51 


0.60 


0.0703 


3840 


8.42 


K5V 


0.67 


0.72 


0.1760 


4410 


7.06 


KOV 


0.79 


0.85 


0.4563 


5150 


5.78 


G5V 


0.92 


0.92 


0.7262 


5560 


5.02 


G2V 


1.00 


1.00 


1.0000 


5790 


4.63 


GOV 


1.05 


1.10 


1.3525 


5940 


4.34 


F5V 


1.4 


1.3 


2.9674 


6650 


3.47 


FOV 


1.6 


1.5 


5.7369 


7300 


2.71 



is generated to have 1000 data points evenly spaced with a 
1-minute cadence centred on the mid-transit time. 

The first noise source we add is shot noise generated by 
taking a random real number from a normal distribution of 
mean zero and standard deviation given by: 

ashot = , , (13) 

where texp is the time between each consecutive measure- 
ment, which is assumed to be one minute in our calculations. 

For the stellar noise, we generate 1000 sinusoidal wave- 
forms with varying periods randomly selected between one 
minute up to 24 hours. These waveforms are then coadded to 
give a single synthetic stellar signal designed to mimic the 
non-stochastic nature of real stellar variability. The r.m.s. 
amplitude of the combined signal is then scaled to 10 ppm, 
which matches the amplitude prediction of BOS for a G2V 
star. 

We consider that the instrument noise is composed of 
hundreds of different noise sources periodic in nature vary- 
ing on timescales from one minute to one day. There is no 
benefit of including timescales longer than this since transit 
events will not last longer than ~ 1 day. The combined r.m.s. 
amplitude of the instrument noise is set to be given by the 
values calculated by Bill Borucki (personal communication) 
and is a function of visual magnitude. For a rriKcp = 12 star, 
we take an r.m.s. instrument noise of 7.47 ppm. 

4.2 Lightcurve fitting 

The noisy lightcurves are generated 10,000 times with the 
correlated noises and photon noise being randomly gener- 
ated in all cases. The lightcurves are then passed onto a 
lightcurve fitting code used to obtain best-fit values for 
T and tc. In all cases we fit for a/R*, i, Rp/R* and tc 
and therefore assume that the out-of-transit baseline is well 
known and its errors are essentially negligible. This is a rea- 
sonable assumption for high quality photometry with large 
amounts of out-of-transit data. 

The fitting code finds the best-fit to the lightcurve 
by utilising a genetic algorithrrQ, PIKAIA (see 



IMetcalfe fc Charbonneaul (|2003l ')'). to get close to a 
minimum in x^- This approximate solution is then used 
as a starting point for a y -minimisatio n performed with 
the AMOEBA routine jPress et all l|l992l )'). The AMOEBA 
solution is tested by randomly perturbing it and refitting in 
20 trials. 

In each subsequent fitting run of the 10,000 lightcurves, 
the AMOEBA routine starts from the original best-fit. The 
resulting values of tc and T are binned and plotted as a 
histogram (figure 2) and t hen compared to the distribution 
expected from the modified lCarter et al.l l|2008l ) expressions. 

4.3 Comparison to the analytic expressions 

For the three cases of a Neptune, a Saturn and a Jupiter, we 
obtain 10,000 estimates of T and tc in each case. We compare 
the distribution of tc and T to the predi cted distribution 
from the expressions of ICarter et al.l l|2008h . In all cases, we 
find excellent agreement between the predicted uncertainties 
and the theoretical values. In figure 2, we plot a histogram 
of the results for the Saturn-case T values, and overlay the 
predicted value of ctt for comparison where the quality of the 
agreement is evident from the plot. Note that this overlaid 
Gaussian is not a fit but a theoretical prediction. 

We find excellent agreement for both or and at^, and 
if anything the theoretical expressions slightly overestimate 
the uncertainties. Thus any results from this study can in- 
fact be considered to be slightly conservative. We prefer to 
adopt a conservative approach in our analysis since there 
may be unexpected sources of noise which increase the tim- 
ing uncertainties. 



5 IMPLEMENTATION OF EXPRESSIONS 

5.1 Earth-transit magnitude limit 

Wo first consider the magnitude limit in detecting an Earth- 
like transit which will allow us to compare the efficiency of 
detecting Earth-like planets and moons. Kepler is designed 
to look at ~ 10^ stars between 6* to 16*'^ mag nitudfEI. For 
each spectral type in table 1, we are able to compute the 
faintest visual magnitude to which a habitable Earth-like 



^ Available from |http: //whitedwarf .org/ parallel/| 



* See |http://kepler.nasa.gov/sci 
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Figure 2. Comparison of the distribution in T found using 
Monte Carlo simulations (scattered points) of a transit of a 
Saturn-like planet versus that from theoretical prediction (smooth 
line). The theoretical prediction for the timing error (42.6s) is 
slightly larger than that found in our simulations (37.3 s). 



transit can be detected. Kepler was conceived with the goal 
of detecting an Earth-Sun transit of 6.5 hours in duration, 
which represents about half the duration if the Earth has 
an impact parameter of zero. However, the maximum mag- 
nitude to which an Earth-like transit could be detected will 
be for the cases of equatorial transits as these maximise the 
transit duration and hence integration time. We establish 
the following criteria for a reliable transit detection: 

• Each transit must be detected to ^ 1-cr confidence. 

• The folded lightcurve must have a significance of ^ S-a. 

• The time between contact points 2 & 3 ^ 1 hour, in 
order to be detected with Kepler's survey-mode cadence of 
30 minutes. 

• At least 3 transits must be detected over the mission 
duration. 

With these criteria, we compute the magnitudes of the 
faintest stars hosting a habitable, transiting Earth which can 
be detected with KCP, for both full transit durations (i.e. 
equatorial transits) and half-durations (i.e. impact param- 
eter chosen to be such that the transit duration is half the 
full transit). The results are given in Table 2. As expected, 
smaller host stars can be fainter due to their smaller ra- 
dius and hence deeper transit depths. We also pick up more 
transits towards smaller stars since the orbital period of the 
habitable Earth decreases. 

Even if Kepler was extended to a mission length of 6 
years, an F5V star would have a habitable zone so far out 
that detecting three transits within this region would be im- 
possible. Hence, we do not consider F type stars in our anal- 
ysis. We also consider G^'-magnitude stars to be the bright 
limit, giving us dmin, the minimum distance for each spec- 
tral type (see Table 2), since the Kepler field is specifically 
chosen to avoid such bright targets. 

5.2 Jupiters vs Saturns vs Neptunes 

We may use the approximate analytic expressions to get 
a handle on the general trends in exomoon detection. The 




Period of planet/years 



Figure 3. Example plot showing the TTV detectability of a 
O.2M9 (^ = 0.4895J exomoon around a Jupiter, Saturn and 
Neptune-like exoplanet for a G2V, m-Kop = 12 star. Values of 
2> 76 are detectable at ^ 8-(j confidence. Saturns present the 
strongest signal due to their low density. The vertical lines repre- 
se nt the stability limit of suc h a moon, calculated using the model 
of \ Barnes & O'BrieA (200 A) . 



first question we may ask is what is the optimum planet to 
search for moons around, out of the three classes of Nep- 
tunes, Saturns and Jupiters? We take each of these planets 
and simulate the detectability of the TTV signal as a func- 
tion of planetary orbital period, Pp, for an O.2M0 exomoon 
with as = 0.4895 Hill radii. We fix the star to be a G2V, 
JTiKop = 12 object and work with i = 90° for simplicity. 

We calculate the signal amplitude and mid-transit time 
uncertainty in all cases and hence find the value for a 
range of orbital periods. In figure 3, we plot as a function 
of host planet's orbital period. The plot reveals that Jupiters 
are the hardest to search for exomoons around whilst Sat- 
urns are the easiest. This is due to Saturn's low density 
meaning a large transit depth but a low enough mass such 
that an exomoon still affects it significantly. We find the 
same order of detectability consistently in many different 
orbital configurations and for TDV as well. 

5.3 ^-rriKcp parameter space 

We now consider an exomoon in the optimum condition of a 
Saturn hosting a single satellite in the habitable zone of the 
host star. The detectability of the timing signals depends 
on the planet-moon separation, as = £,Rh, the apparent 
magnitude of the host star in the Kepler bandpass, m-Kop, 
mass of the exomoon, Ms, and finally the host star's mass 
and radius (i.e. spectral type). Note that Rh denotes the 
Hill radiufl 

TDV increases with lower values of ^ and TTV increases 
with higher values of ^, but we must maintain a balance 
such that at least one of the effects is detected to 8-(j and 
the other to 3-(t confidence (dictated by our previous de- 
tection criteria). To find the limits of interest, we assume a 
zero-impact-parameter transit of our optimum planet-moon 

^ Not e that we denote the quantity asJRH as as opposed to x 
in [Kipping ( 2009i) 

and iKip ping 12009bl 'l to avoid confusion with 
the chi square distribution. 
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Table 2. Faintest stars for which a habitable transiting Earth could be detected for different star types. Final column gives maximum 
distance of such a star with a full transit duration (Df^u), based on absolute magnitude in the Kepler bandpass. A blank indicates that 
no magnitude can satisfy the detection criteria. 



Star type 


Dhalf "iKcp.max 


D full ™'Kep,max 




^max 1 


M5V 




18.111 


0.68 


179.56 


M2V 


15.992 


16.190 


2.00 


218.78 


MOV 


15.206 


15.465 


3.28 


256.45 


K5V 


14.278 


14.618 


6.14 


324.79 


KOV 


13.286 


13.712 


11.07 


385.83 


G5V 


12.762 


13.240 


15.70 


440.56 


G2V 


12.236 


12.763 


18.79 


423.25 


GOV 


11.520 


12.106 


21.48 


357.44 


F5V 










FOV 
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be ~ 0.3 Ganymede masses, assuming the maximum pro- 
grade orbital distance is 0.4895 Hill radii, as calculated by 
iDomingos et al.l (|2006l ). In contrast, an M2V star allows for 
a habitable Saturn to hold onto a 0.4-Earth-mass moon for 
over 5 Gyr. 

A second lower limit exists from tidal forces and the 
Roche limit. An upper limit is given by the fact we require 
three transits in a 4 year observation duration and thus the 
most distant habitable zone assumed here corresponds to a 
period of 1.33 years, which excludes the FOV and F5V stellar 
types. 



Figure 4. Detectable range of a habitable -^M^ exomoon for an 
MOV star with respect to apparent magnitude and ^ parameter 
space, with Heplev-class photometry. The grey area represents 
the area satisfying our detection criteria of detecting TTV or 
TDV to 8-cr and the other timing effect to 3-cr, shown by the 
four curves. Additional constraints are the lower Roche limit, 
the upper dynamical stability limit and the bright-star magnitude 
limit. 



arrangement. The orbital period of the host planet is always 
given by Phab- We prefer to calculate this best-case scenario 
since we may then consider it unfeasible to detect habitable 
moons beyond this limit. 

For a given star type and exomoon mass, we plot the 
contours of rriKep and ^ which provide TTVs of 3- and 8- 
o" confidence and then repeat for TDVs of the same confi- 
dences. There are two possible acceptance criteria: i) TTV 
confidence is ^ 3-cr and TDV confidence is ^ 8-a; ii) TTV 
confidence is ^ 8-a and TDV confidence is ^ 3-a. The loci 
of points below these lines represents the potential detec- 
tion parameter space. For each star type and exomoon mass, 
these loci will be different. In figure 4 we show a typical ex- 
ample for an MOV- type star and Ms ~ ^M^. 

There are certain cases wh i ch con strain the stars of in- 
terest. The I Barnes fc O 'BrienI (|2002l ) limit may be calcu- 
lated for a Saturn harbouring an exomoon in a system of 5 
Gyr age. This suggests that the maximum moon mass that a 
habitable zone Saturn could hold around a M5V star would 



5.4 Magnitude quartiles 

For figure 4, we are able to convert the plot into three key 
numbers by assuming that the probability distribution of 
exomoons with respect to ^ is approximately flat. This sim- 
ple a-priori approximation allows us get a handle on some 
typical magnitude limits. This assumption effectively con- 
verts the two-dimensional array of ^i-mKGp,i into a one- 
dimensional list of just rriKop.i. We then simply take the 
quartiles of this list to obtain estimates for the 25%, 50% 
and 75% catch-rate values of rriKcp. Physically speaking, we 
are saying that, for example, at the 50% catch-rate value of 
rriKcp, there is a 50% probability of detecting an exomoon 
if an exomoon exists with a flat prior distribution between 
^min and ^ 

max . 

Accordingly, we can calculate the magnitude limits to 
detect 25%, 50% or 75% of the exomoons in the given sam- 
ple, which we label as the upper, middle and lower quartiles 
respectively. For the case shown in figure 4, the quartile val- 
ues are mxep ~ 12.0, 11.5 and 10.9 respectively. 

We may then calculate these quartile values for different 
exomoon masses, for a given star type. Effectively, we are 
able to determine the minimum detectable exomoon mass as 
a function of magnitude for three different detection yields: 
i) 25%, ii) 50%, and iii) 75%. In figure 5, we show an example 
of one of these plots for an M2V star with a habitable-zone 
Saturn host planet. We do not consider magnitudes brighter 
than 6"^ magnitude, since Kepler's field-of-view has been 
selected to omit such stars. 
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Figure 5. Minimum detectable exomoon mass around a Saturn- 
like planet host in the habitable zone of an MOV star. The top 
line is the 75% catch-rate, the middle is 50% and the lowest is for 
a 25% exomoon catch-rate. We also mark three key magnitudes 
of interest; 9, 12 and I4 magnitudes. 




Mass of host star/M^ 



Figure 6. Minimum detectable exomoon mass around a Saturn- 
like planet host in the habitable zone of various star types. Con- 
tours show the different magnitude limits for a 25% exomoon- 
catch-rate a nd overlaid is the mass stability limit (black, solid), 
found usina ^Bames & O'Brien ij200k) . 



5.5 Minimum detectable habitable exomoon 
masses 

To complete the picture, we wish to see how these values dif- 
fer for various star types. We repeat the analysis performed 
above for M5V, M2V, MOV, K5V, KOV, G5V, G2V and 
GOV-type stars. In figure 6, we plot the minimum detectable 
exomoon mass, for several contours of visual magnitude, as 
a function of stellar mass, in the 25% detection-yield case. 
We also show the upper l imit on moon mass calculated from 
iBarnes fc O'BrienI ([20o3 ) but this is only an approximation 
since we have assumed values of Qp and k2p. 

The mKep = 12.5 limit almost intersects the 

iBarnes fc O'BrienI (|2002l 'l limit at O.4M0 and ~ 0.4M(b for 
the exomoon corresponding to an M2V star ~ 40 pc away. 
For comparison, within 10,000pc'' (a spherical shell of ra- 



Figure 7. Distance limit for detecting a 1 habitable exomoon 
(solid) and a 1 exoplanet (dot-dashed) as a function of stellar 
mass, for main sequence stars. We assume a Saturn-like planet 
as the host for the exomoon case. We also mark the lower stellar 
mass stability limit of ~ 0.475 Mq and the upper mass limit of 
~ 1 Mq imposed by the requirement to observe ^ 3 transits. 



dius 13.4 pc) there are ~ 800 ma i n-sequ ence stars, ~ 630 of 
which are M-dwarfs (see lLedrewl l|200ll )'). 

If we set a lower limit of 10 pc for our target star, then 
the lowest-mass habitable-zone exomoon which could be de- 
tected with KCP around a M5V star would be 0.09M©, 
which is above the maximum allowed stable mass limit of 
O.OOSMq. However, moving up to an M2V star we find a 
minimum detectable habitable exomoon mass of O.ISM^ 
which is below the stability limit of O.41M0. We therefore 
conclude that the minimum detectable habitable-zone exo- 
moon mass with KCP is ~ 0.2Mq. 



5.6 1 M© Limit 

For the 25% detection yield case, we can convert the mag- 
nitude limit for detecting a 1 habitable exomoon into a 
distance limit by making use of the absolute magnitudes for 
each star type. The distance limit for detecting a 1 hab- 
itable exomoon may then be compared to the limit found for 
a transiting 1 Rq habitable exoplanet. In figure 7, we com- 
pare the two distance limits, which reveals that the distance 
limit for moons takes the same shape as the transit limit 
but is ~ 1.6 times lower. This translates to a volume space 
diminished by a factor of ~ 4. 

It is not our intention to estimate accurately how many 
observable stars might be contained within this volume, but 
one can deduce an order-of-magnitude estimate. Assuming 
that there are about 10^ useable stars in Kepler's field of 
view, we estimate roughly that 25,000 stars would be within 
range for habitable-exomoon detection. Extrapolating Ke- 
pler's field of view to a band centred on the Galactic plane 
between b = ±15°, and simplistically assuming the same 
number of target stars per square degree, we obtain about 
2 X 10^ stars that are accessible to KCP for a search for 
habitable-zone exomoons, if each star had a habitable-zone 
gas giant similar to Saturn. 
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5.7 Detectability of an Earth-Moon analogue 

Finally, we consider the detectability of an Earth-Moon sys- 
tem. Using the orbital parameters of the Sun-Earth-Moon 
system and assuming the most favourable configuration of 
i — 90°, we find the Earth would exhibit a TTV r.m.s. am- 
plitude of 112.4 seconds and a TDV of 13.7 seconds. The 
TDV is particularly low due to the large value of ^ for the 
Moon. 

In comparison, KCP of a y = 6 star would yield timing 
errors of at^ — 389.6 seconds and ctt = 779.2 seconds. Even 
for the most favourable case then, the TTV has a signal-to- 
noise of 0.3, which would be quite undetectable. We there- 
fore conclude that KCP cannot detect a Sun-Earth-Moon 
analogue through transit timing effects. 



6 CONCLUSIONS 

Our analysis finds that habitable exomoons are detectable 
up to ~ 100-200 pc away around early-M, K and later-G 
dwarf stars with the Kepler Mission or photometry of equal 
quality (KCP). This photometric quality should be sensitive 
down to 0.2M(^ habitable exomoons in the idealised cases 
and could survey ~ 10^ stars for IM® habitable exomoons, 
with a full Galactic Plane survey. This number is around 
25,000 stars for Kepler's field-of-view. 

We find that Saturn-like planets are the ideal host can- 
didates for detection due to their large radius to mass ratio. 
Additionally, we find lower-mass exomoons may be found 
around M-dwarfs due to the closer habitable zone permit- 
ting a larger number of transits in a 4-year window. 

The exciting prospect of discovering a habitable exo- 
moon is well within the grasp of KCP and whether such 
worlds exist or whether they will be classed as truly hab- 
itable worlds are questions we can hope to answer in the 
coming years. It is interesting to note that ~ O.2M0 is the 
minimum habitable mass for a planet or m oon used by sev- 
eral ot her sources in the literature, in cluding Ra ymond et al.l 
(|2007D (0.3M®) and lWilhams et al I (l997) (0.12Mg3). It has 
also been proposed bv lScharj l|2006l ) that the habitable zone 
could be extended for exomoons due to tidal heating and so 
our calculations may infact be an underestimate. 

We find that Kepler will be incapable of find- 
ing Ganymede-mass moons around Saturn-like planets. 
ICanup fc WardI l|2006l ) have suggested that the maximum 
moon mass which could form from a planetary debris disk is 
2 X 10"* of that of the primary's and therefore, if this rule 
holds true, it unlikely KCP would detect any moons which 
formed around a planet in such a way. However, moons 
ab ove this mass limit are st ill dynamically stable (according 
to iBarnes fc O'BrienI l|2002h ) and could have been captured 
by a planet or formed through an impact, for example like 
the formation of Triton and the Moon respectively. Whether 
or not such objects are common is unknown but KCP could 
make the first in-depth search. 

Our results suggest it is easier to detect an Earth-like 
exoplanet than an Earth-like exomoon around a gas giant. 
However, we have no statistics to draw upon to estimate 
which of these scenarios is more common. If a roughly equal 
number of both are discovered, it would indicate that the 
latter is more common due to the detection bias. 



These results highlight the promising opportunity of 
making the first exomoon detection using the Kepler tele- 
scope, or photometry of equivalent quality, especially the 
feasibility of detecting habitable-zone exomoons. All-sky 
surveys focussing on bright M-dwarf stars would be ideally 
placed to search for habitable exomoons in greater depth 
and thus a telescope like TESS could continue the search 
after the Kepler Mission ends. 
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